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Abstract. We give an overview of the current status of investigations of the polarization of 
gluons in the nucleon. We describe some of the physics of the spin-dependent gluon parton dis- 
tribution and its phenomenology in high-energy polarized hadronic scattering. We also review 
the recent experimental results. 

1. Introduction 

For many years now, spin has played a very prominent role in QCD. The field of QCD spin 
physics has been driven by the hugely successful experimental program of polarized deeply- 
inelastic lepton-nucleon scattering (DIS) [1]. One of the most important results of this program 
has been the finding that the quark and anti-quark spins (summed over all flavors) provide only 
about a quarter of the nucleon's spin, AS ~ 0.25 in the proton helicity sum rule [2j [3J 0] 

\ = ^AS(Q 2 ) + AG(Q 2 ) + L q (Q 2 ) + L g (Q 2 ) , (1) 

implying that sizable contributions to the nucleon spin should come from the gluon spin 
contribution AG(Q 2 ), or from orbital angular momenta L qjg (Q 2 ) of partons. Here, Q is the 
resolution scale at which one probes the nucleon. The Q 2 -dependence of the various contributions 
to the proton spin is predicted in perturbative QCD through its evolution equations [H EJ E]. 
To lowest order (LO), the quark and anti-quark spin contribution AS(Q 2 )/2 does not depend 
on Q 2 . Figure [T] shows a LO toy calculation of the Q 2 -evolution of the contributions to the 
proton spin in Eq. ([1]), assuming that at an initial scale Qq = 1 GeV we have AS = 0.25, 
AG = L q = 0.2, L g = -0.025. The rise of AG oc log(Q 2 ) or l/a s (Q 2 ) (compensated by an 
opposite evolution of L g ) is an important prediction of QCD and awaits experimental testing. 
For the initial conditions chosen here, the evolution leads to large positive values of AG. We 
note that at asymptotic Q 2 the total quark and gluon angular momenta, ^AT, + L q and AG+L g , 
respectively, become roughly equal [I]. 

To determine the gluon spin contribution on the right-hand-side of Eq. ([T]) has become a major 
focus of the field. Like AS, it can be probed in polarized high-energy scattering. Several current 
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Figure 1. LO toy calculation of the Q 2 - 
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experiments are dedicated to a direct determination of the spin-dependent gluon distribution 
Ag(x,Q 2 ), 

Ag(x,Q 2 ) =g + (x,Q 2 )-g-(x,Q 2 ) , (2) 

where g + (g~) denotes the number density of gluons in a longitudinally polarized proton with 
same (opposite) sign of helicity as the proton's, and where x is the gluon's light-cone momentum 
fraction. The field-theoretic definition of Ag is 



written in A + = gauge. G^ v is the QCD field strength tensor, and G^ u its dual. The integral 
of Ag(x, Q 2 ) over all momentum fractions x becomes a local operator only in A + = gauge and 
then coincides with AG(Q 2 ) [21 Ej. The COMPASS experiment at CERN and the HERMES 
experiment at DESY attempt to access Ag(x,Q 2 ) in charm- or high-p-j- hadron final states in 
photon-gluon fusion j*g — > qq. A new milestone has been reached with the advent of the first 
polarized proton-proton collider, RHIC at BNL [8, 9j. RHIC will provide precise and detailed 
information on Ag, over a wide range of x and Q 2 , and from a variety of probes. 

2. Model estimates of Ag 

Before we discuss in some detail the phenomenology of Ag in polarized high-energy scattering, let 
us briefly address some of the available theoretical expectations for Ag and its integral. As was 
first pointed out in [10], it is possible to estimate the operator matrix element corresponding 
to AG in non-relativistic quark and bag models. In such models, for example, baryon mass 
splittings result from lowest-order exchange of transverse gluons, and the associated forces 
are spin-dependent. One obtains estimates [TTJ Q2] for AG(Q 2 « 1 GeV 2 ) of about 0.2 to 
0.3. In a sense, these are "natural" values since they are of the order of the proton spin 
itself. Very recently, for the first time model calculations of the x-dependence of Ag have 
been presented |12j . The resulting distribution is positive everywhere and of moderate size. 
The more and more precise experimental constraints on Ag will likely motivate further model 
investigations, which ultimately might lead to new insights into QCD. Likewise, it is to be hoped 
that lattice calculations, which are becoming ever more powerful, will be able to address gluonic 
observables in nucleon structure in the future [13]. 

Other considerations, based in part on perturbation theory, led to the prediction of a very 
large gluon polarization in the nucleon. The peculiar evolution pattern of AG(Q 2 ) oc l/a s (Q 2 ) 




Aff(x,Q 2 ) = — / d\e lXxP+ (P,S\G + »(0)G + v (\n)\P,S) 



visible in Fig. [T] inspired ideas [13] that a reason for the experimentally found small size of the 
proton's axial charge should be sought in a "shielding" of the quark spins due to a particular 
perturbative part of the DIS process j*g — > qq. The associated cross section is of order 
a s (Q 2 ), but the (Revolution of AG(Q 2 ) would compensate this suppression. We note that this 
interpretation of the axial charge, however, corresponds to a particular choice of factorization 
scheme. To be of any phenomenological relevance, such "anomalous" models would require a 
very large positive gluon spin contribution, AG > 1.5, even at a low scale of 1 GeV or so. As 
we shall see below, initial experimental data now appear to make such a scenario very unlikely. 



3. Ag and scaling violations in polarized DIS 

In principle, a clean determination of Ag(x, Q 2 ) is possible by investigating scaling violations of 
the spin-dependent proton structure function g\(x,Q 2 ) which is measured in polarized DIS. To 
leading order of QCD, g\ can be written as 



gi{x,Q 2 ) = \Y,e 2 q [Aq{x,Q 2 ) + Aq{x,Q 2 )] , (4) 



where the Aq and Aq are the quark and anti-quark helicity distributions. QCD predicts the 
Q 2 -dependence of the densities through the spin-dependent DGLAP evolution equations [5]: 

d ( Ag \ 2 _ f 1 dz f AP qq (a s (Q 2 ),z) AP qg (a s (Q 2 ),z)\ ( Aq \ fx , 
dhxQ 2 \Ag)^ } ~ J x z \AP gq {a s {Q 2 ),z) AP gg (a s (Q 2 ), z) )\ Ag ) \ z ' y 

(5) 

the APij are the spin-dependent "splitting functions" [H [15] which are evaluated in QCD 
perturbation theory. As one can see, Ag contributes to the scaling violations of g\. Nonetheless, 
Ag has been left virtually unconstrained (see, for example, [TBJ [T71 [TBI [H2 [20]) by the scaling 
violations observed experimentally in polarized DIS. This is due to the very limited lever arm 
in Q 2 of the fixed-target experiments. Figure [2J shows current theoretical "uncertainty bands" 
for Ag from DIS scaling violations. At best, a tendency toward a positive Ag is seen. We 
note that a recent new analysis by the COMPASS collaboration [2T] using their latest deuteron 
DIS data [22J finds two "allowed" regions for Ag, one with positive, one with negative gluon 
polarization. Clean and precise extractions of Ag(x,Q 2 ) over a wide range of x and Q 2 from 
scaling violations of g\ would become possible at a polarized electron-ion collider, EIC [23], 
thanks to its vastly larger kinematic reach. 
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t igure z. Results for xAg(x, Q 2 = 5 GeV 2 ) 
from several analyses [HI [T71 [18] of polarized 
DIS. The various bands indicate ranges in Ag 
that were deemed consistent with the DIS 
scaling violations in these analyses. From [9]. 



4. Access to Ag in polarized proton-proton scattering at RHIC 

The measurement of gluon polarization in the proton is a major focus and strength of RHIC [HI [9]. 
The basic concept that underlies most of spin physics at RHIC is the factorization theorem [23] . 
It states that large momentum-transfer reactions may be factorized into long-distance pieces that 
contain the desired information on the spin structure of the nucleon in terms of its universal 
parton densities, and parts that are short-distance and describe the hard interactions of the 
partons. The latter can be evaluated using perturbative QCD. As an example, we consider the 
double-spin asymmetry for the reaction pp — ► irX, 

a ++ - a + - Act , , 

a ++ + a + a 

where the superscripts denote helicities of the initial protons. We assume the pion to be produced 
at high transverse momentum pr, ensuring large momentum transfer. Then, up to corrections 
suppressed by inverse powers of pr- 

dAo = J2Af a ® Af b ® dAa c ab Dl (7) 

abc 

for the polarized cross section, where (g> denotes a convolution. The A/j are the polarized 
parton distributions, and the pion fragmentation functions. The sum in Eq. (JJJ) is over all 
contributing partonic channels a+b — * c+X, with dAa^ b the associated spin-dependent partonic 
cross section. In general, a leading-order estimate of ([7j) merely captures the main features, but 
does not usually provide a quantitative understanding. Only with knowledge of the next-to- 
leading order (NLO) QCD corrections to the dAa^ can one reliably extract information on the 
parton distribution functions from the reaction. 

Several different processes will be investigated at RHIC [HI [9] that are very sensitive to 
gluon polarization: high-p^ prompt photons pp — ► 7X, jet or hadron production pp — ► jet A, 
pp — ► hX, and heavy-flavor production pp — > (QQ)X. An important role for the determination 
of Ag will be played by measurements of two-particle, jet-jet (or hadron-hadron) and photon- 
jet correlations. For these, at the leading order approximation, the hard-scattering subprocess 
kinematics can be calculated directly on an event-by-event basis, giving an estimate of the gluon 
momentum fraction [25]. In addition, besides the current y/s = 200 GeV, also y/s = 500 GeV 
will be available at RHIC at a later stage. All this will allow to determine Ag(x, Q 2 ) in various 
regions of x, and at different scales. Essentially all tools are in place now for treating the 
spin-dependent reactions relevant at RHIC at NLO [251 1271 1251 129] . 

We emphasize that there have already been results from RHIC that demonstrate that the NLO 
framework is very successful. Figure [3] shows comparisons of data from PHENIX and STAR for 
single-inclusive cross sections for ir° [30] [31] , jets [32] and photons [33J with corresponding NLO 
calculations [26] [28] [29] [53]. As can be seen, the agreement is overall excellent. We note that 
an agreement between data and NLO calculations like the one seen in Fig. [5Jis not found in the 
fixed-target regime [35] (it has recently been shown that in this regime large logarithmic terms 
at yet higher orders are important and need to be resummed for a more successful theoretical 
description [36] ). In Fig. [3J we decompose the NLO mid-rapidity tt° cross section into the 
relative contributions from the various two-parton initial states [9]. It is evident that processes 
with initial gluons dominate. 

The results shown in Fig. [5] give confidence that the theoretical NLO framework may be used 
to determine the spin-dependent gluon density from RHIC data. Results for All in PP — > ttA are 
now available from PHENIX [37J, and All for single-inclusive jet production has been measured 
by STAR [52]. The results are shown in Fig. [5] The curves shown in Fig. [5] represent the All 
values calculated at NLO for a range of gluon distributions from [T7], from a suggested very large 
positive gluon polarization ("GRSV-max") with an integral AG = 1.9 at scale Q = 1 GeV, to a 
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Figure 3. Data for the cross section for single-inclusive ir production pp — > tt X at 
^fs = 200 GeV at mid-rapidity from PHENIX (upper left, [30]) and at forward rapidities from 
STAR (lower left, [31]), for mid-rapidity jet production from STAR (upper right, [32J), and for 
mid-rapidity prompt-photon production from PHENIX (lower right, |33j). The lines show the 
results of the corresponding next-to-leading order calculations [261 [28| [29]. 
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"maximally" negative gluon polarization, ( u Ag = —5"), for which AG(lGeV 2 ) = —1.8. These 
two distributions span the "GRSV-band" shown in Fig. [2l The curves labeled "GRSV-std" 
represent the best fit of [T7] to the polarized DIS data (solid line in Fig. [2]), which has a more 
"natural" AG(1 GeV 2 ) of about 0.4, and the results for 11 Ag = 0" correspond to very little gluon 
polarization, AG(1 GeV 2 ) = 0.1. One can see that the data are already discriminating between 
the various Ag distributions. A very large gluon distribution, as proposed in the context of 
the "anomaly scenario" (see discussion above) and corresponding roughly to the curves labeled 
"GRSV max" , appears to be strongly disfavored. 
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Figure 5. Data for the double-spin asymmetry for mid-rapidity single-inclusive n production 
at y/ii = 200 GeV from PHENIX [37] (left), and for jet production from STAR [32] (right), 
compared to NLO predictions for several polarized gluon distributions of [17j . 

Figure [6] shows NLO predictions [28] for the double-longitudinal spin asymmetry All for the 
reaction pp — > 7V at RHIC, based on the "gluon uncertainty" band displayed in Fig. [2j Prompt 
photons are much less copiously produced than pions at RHIC, resulting in larger statistical 
uncertainties. The measurement of this asymmetry will therefore take some time at RHIC. 
Nonetheless, the reaction pp — ► 7X is of great importance because of its direct sensitivity to Ag 
through the clean "Compton-like" process qg — ► 79. The spin asymmetry for this reaction is 
linear in Ag and therefore directly determines the sign of the distribution. The plot also shows 
the experimental uncertainties expected at RHIC (PHENIX) for 65/pb collected luminosity [9]. 
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Figure 6. "NLO theory band" for 
mid-rapidity prompt-photon production at 
PHENIX. The band illustrates the current 
"uncertainty" due to Ag, using the GRSV Ag 
band shown in Fig. El The "errors" are pro- 
jections. Taken from [9]. 



5. Ag from photon-gluon fusion 

A way to access Ag in lepton-nucleon scattering is to measure final states that select the photon- 
gluon fusion process. These are heavy-flavor production, Ip — ► ccX, and single- or di-hadron 
production, Ip — ► hX or £p — > h\h<iX, where the hadrons have large transverse momentum. 
Figure [7] compiles the current results [2 H 138 1 f39l I40j for extractions of Ag from these reactions. 
We note that, unlike at RHIC, the success of the perturbative-QCD hard-scattering description 
has not been established for these observables in the kinematic regimes of interest here. Also, 
the translation of the measured spin asymmetry into Ag at a certain single momentum fraction, 
currently only possible at leading order, is fraught with large uncertainties. 
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Figure 7. SMC [38], COMPASS [39], and 
HERMES [30] results for gluon polarization, 
extracted from photon-gluon fusion. Taken 
from [21J. 



6. Global Analysis 

The eventual determination of gluon polarization will require consideration of all existing data 
through a "global analysis" that makes simultaneous use of results for all probes, from RHIC 
and from lepton scattering. The technique is to optimize the agreement between measured spin 
asymmetries, relative to the accuracy of the data, and the theoretical spin asymmetries, by 
minimizing the associated x 2 function through variation of the shapes of the polarized parton 
distributions. The advantages of such a full-fledged global analysis program are manifold: (1) 
The information from the various reaction channels is all combined into a single result for Ag(x). 
(2) The global analysis effectively deconvolutes the experimental information, which in its raw 
form is smeared over the fractional gluon momentum x, and fixes the gluon distribution at 
definite values of x. Figure [8] highlights the importance of this. The figure shows [H] the 
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Figure 8. NLO dAa / dpxdlog^x) (arbitrary normalization) for the reaction pp — » tt°X at 
RHIC, for p T = 2.5 GeV and six different values for AG(p 2 ) at p » 0.4 GeV [17]. The shaded 
areas denote in each case the re-range dominantly contributing to dAa. From [41j . 



contributions of the various regions in gluon momentum fraction to the mid-rapidity spin- 
dependent cross section for pp — > ir°X at RHIC, for six different sets of polarized parton 
distributions [17] mostly differing in the gluon distribution. The pion's transverse momentum 
was chosen to be 2.5 GeV. One can see that the distributions are very broad, and that the 
x-region that is mostly probed depends itself on the size and form of the polarized gluon 
distribution. This makes it very difficult to assign a good estimate of the gluon momentum 
fraction to a data point at a given pion transverse momentum. The global analysis solves this 
problem. 

The further advantages of a global analysis are: (3) State-of-the-art (NLO) theoretical 
calculations can be used without approximations. (4) It provides a framework to determine 
an error on the gluon polarization. (5) Correlations with other experiments, to be included in 
X 2 and sensitive to degrees of freedom different from Ag, are automatically respected. Global 
analyses of this type have been developed very successfully over many years for unpolarized 
parton densities. Examples of early work on global analyses of RHIC-Spin and polarized DIS 
data in terms of polarized parton distributions are |42 } 143 } l4"4"]. 



7. Conclusions and Outlook 

While the initial data from RHIC and from the dedicated studies in lepton scattering shown 
above point to a small or moderate size of the gluon polarization in the x-region currently 
accessible, statements about the gluon contribution to the proton spin, AG, are really not 



possible yet and will require the global analysis just described. A crucial issue will eventually 
be the behavior of the extracted Ag(x) at the smallest accessible x, which are reachable in 
500 GeV running at RHIC and in correlation studies involving final states produced at forward 
angles |25j . It is possible that a significant contribution to AG comes from relatively small x. 
As one example, we show in Fig. [9] the "running integral" J x . dxAg(x,Q 2 ) at Q 2 = 10 GeV 2 , 
normalized to the full integral AG(Q 2 ), for the gluon distribution in the NLO GRSV "standard" 
set |17j . As one can see, at this scale about 30% of the integral come from x < 10~ 2 . RHIC 
will likely be able to constrain A g down to values somewhat smaller than that, but the example 
shows that it might eventually be necessary to push to x < 10~ 3 and below. This could be 
achieved at a high-energy polarized electron-proton collider [23J. 
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